Guiding chemical pulses through geometry: Y-junctions 
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We study computationally and experimentally the propagation of chemical pulses in complex 
geometries. The reaction of interest, CO oxidation, takes place on single crystal Pt(llO) surfaces 
that are microlithographically patterned; they are also addressable through a focused laser beam, 
manipulated through galvanometer mirrors, capable of locally altering the crystal temperature and 
thus affecting pulse propagation. We focus on sudden changes in the domain shape (corners in a 
Y-j unction geometry) that can affect the pulse dynamics; we also show how brief, localized temper- 
ature perturbations can be used to control reactive pulse propagation. The computational results 
are corroborated through experimental studies in which the pulses are visualized using Reflection 
Anisotropy Microscopy. 

PACS numbers: 82.40.Ck, 82.40.Np, 05.45.-a 



I. INTRODUCTION 

Reaction control is an essential task in the chemical 
process industry; it can stabilize unstable (but profitable) 
operating states, and, in a multiobjective context, can 
enhance reactor performance while satisfying safety and 
environmental constraints. Non-steady state (e.g. peri- 
odic) operation may result in better average productivity 
or selectivity than steady state operation (as suggested 
in the early 1960s in the pioneering work of F. Horn (see 
Ref. Q and references therein)). At a much finer scale, 
the nonlinear characteristics of heterogeneous catalytic 
reactions on single-crystal catalysts could be exploited to 
improve overall yield and selectivity through spatiotem- 
porally resolved actuation 0, [3| . Micro-lithography, by 
designing the shape of reactive domains as well as pre- 
scribing the geometry and statistics of heterogeneous in- 
clusions, provides a different avenue of "talking" to the 
local dynamics of catalytic reactions (e.g. Ref. 0|). The 
phenomenology of pattern formation on micro-composite 
catalytic surfaces can be much richer than that in uni- 
form media. There have been extensive experimental and 
computational studies in this area, both from our group 
11, lijl EJi^ [H [H and in the work of oth- 
ers [13, [1^ [13 • For the low pressure CO oxidation on 
Ft (110), the heterogeneities deposited onto the catalytic 
Ft surface can be either inert, such as Ti02, or may 
consist of different active catalysts for the reaction, e.g. 
metals like Rh and Fd. Such fields of heterogeneous in- 
clusions affect the reaction dynamics and the formation 
of patterns on the catalytic Ft surface largely through 
their interfaces. 

The heterogeneity in such composite catalysts is solely 
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spatial; their geometry is determined upon construction 
and does not change with time. In more recent stud- 
ies, spatiotemporal forcing (using a galvanometer mirror- 
manipulated focused laser beam) was applied during CO 
oxidation on a Ft (110) single-crystal [3. Fulses and 
fronts, the basic building blocks of patterns, could be 
formed, guided and destroyed with an addressable laser 
beam; the addressability function was also used to en- 
hance catalytic reactivity in Ref. 3]. 

From the implementation of chemical "logical gates" 
(e.g. Ref. 18]) to the recent spatiotemporal control of 
morphogenetic patterns in drosophila embrvos |19l. and 
from drop formation control in microfluidics [2^ to the 
guidance of matter waves in Bose-Einstein condensates 
|21j . guiding pulses and fronts in complex geometries 
is an essential task in spatiotemporal pattern control. 
What makes it increasingly possible is the combination 
of spatiotemporally finely resolved sensing combined with 
spatiotemporally finely resolved actuation. Whether the 
pattern- forming system is naturally photosensitive, like 
the chemical waves in Ref. or rendered photosensi- 
tive, like the neurons in Ref. through genetic ma- 
nipulation, rapid technological developments in finel y re - 
solved optical actuation techniques (e.g. see Ref. [2^ ) 
are rapidly and radically changing the experimental ex- 
ploration and control of spatiotemporal pattern forma- 
tion. 

In this paper we explore the combined effects of ge- 
ometry and of spatiotemporally resolved addressability. 
In particular, we explore the influence of inert bound- 
ary geometry on the dynamics of propagating pulses for 
CO oxidation on a micro-composite Ft (110) surface, and 
show how these dynamics can be altered using an ad- 
dressable laser beam. Our geometry of choice for pulse 
propagation studies is a Y-shaped junction structure, in 
which sharp boundary curvature changes (corners) can 
effectively dictate pulse propagation. We systematically 
explore the effect of varying the junction geometry on 



reactive pulse propagation; a related rhomb geometry is 
also studied with qualitatively similar results. We then 
show how one can actively alter the phenomena dictated 
by geometry through the use of single-shot spatiotempo- 
rally localized laser heating. The computational predic- 
tions are validated by experimental observations of reac- 
tive pulses visualized through Reflection Anisotropy Mi- 
croscopy (RAM) 25J. 

II. MODELING 

In our simulations we use the three- variable Krischer- 
Eiswirth-Ertl reaction-diffusion model for CO oxidation 
on Pt(llQ) with a surface phase transition described 
in Ref. [2^. This surface reaction follows a Langmuir- 
Hinshelwood mechanism: 

CO + * ^ COads 
2 * +O2 ^ 20ads 
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FIG. 1: (Color online) Pulse propagation in different Y- 
j unction structures. By choosing appropriate geometric pa- 
rameters, we are able to dictate the pulse transmission pat- 
terns, (a) /i = 5, ^ = TT, (b) h — 3.7, ^ = tt, (q) h — 3.5, 
^ = TT, (d) /i = 3.5, d — Itt. Other parameters: = 5, 
w 1.93, a = ^TT, T = 535. 5i^, Pco 4.95 x 10"^m6ar, 
P02 = 2.0 X \^-^mhar. 



COads + Oads ^ 2 * +CO2 T 

accompanied byalx2^1xl phase transition of 
the Pt(llO) surface due to CO adsorption. When the 
coverage of CO lies between 0.2 and 0.5, the fraction of 
1x1 surface increases monotonically as the CO coverage 
increases. The sticking coefficient of oxygen is 50% higher 
on the 1x1 surface as compared to the 1x2 surface. 
The 1x1 phase favors oxygen adsorption, which leads to 
reactive consumption of CO. This can lead to oscillatory 
behavior of the reaction, and also allows the formation 
of propagating reaction pulses. 

The equations for this kinetic model are 




klU — k2uv + DuV u 
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where hy v and w we denote the surface coverage of 
CO and O, and the surface fraction of 1 x 1 phase re- 
spectively. The adsorption rate constants for CO and 
O2, ku and ky respectively, are considered to be constant 
within the temperature range considered here. The rate 
constants /ci,/c2 and ks for the desorption, reaction and 
surface phase transition are given by the Arrhenius for- 
mula ki = {k^)iexp{—Ei/RT) and T is the temperature 
of the single-crystal. The function f{u) is fit to experi- 
mental data to give the rate of surface phase transition 
as a function of li, the CO surface coverage, as follows: 
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0.2 <u< 0.5 
u > 0.5 



For simplicity, the diffusion of CO is taken to be 
isotropic and no-fiux boundary conditions are used in 
our simulations. To reflect the influence of laser heating 
on the dynamics of the reaction, we approximate the lo- 
cal temperature increase caused by a laser spot through 
a local Gaussian temperature field. The heat generated 
by the reaction (~ ImW) can be neglected compared to 
the power of the laser beam(~ IW). Since the diffusiv- 
ity of adsorbed CO (~ 10~^cm^/s) is much smaller than 
the thermal diffusivity constant of Pt (~ 10~ ^cm^/s), 
we can assume the local Gaussian temperature field is 
established (resp. vanishes) instantaneously as the laser 
beam is applied (resp. removed) 23,(2^. In our simu- 
lations, we use the commercial Finite Element package 
FEMLAB to compute the time-dependent evolution of 
reactive pulses in 2D Y-j unction geometries; our meshes 
typically consisted of ~ 12,000 linear elements. 



III. COMPUTATIONAL RESULTS 

A. Pulse propagation in a Y-junction 

The geometry of the Y-junction structure is described 
through the parameters VK, /i, a, and 6 (see the first 
snapshot in Fig. da)). When a reactive pulse reaches 
the corners of the Y-junction that are denoted by the 
two angles a and ^, the pulse front starts becoming con- 
vex. The local travelling speed of the pulse decreases due 
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FIG. 2: Geometry of a corner structure, r = 



to this convex curvature . When the local curva- 

ture required to "go around" one of the corners exceeds a 
certain limit, the pulse loses stability, "decollates" from 
the boundary, and disappears |31^ . With appropriate an- 
gles a and ^, by adjusting the position of the "prow" 
to the exit (i.e. adjusting /i), we can force the pulse to 
choose none, one, or both of the channels to propagate in 
(Fig.Q). The behavior of the pulse also depends strongly 
on the selected partial pressures of CO and O2. Under 
reaction conditions for Fig. Q decreasing the amount of 
CO in the gas phase stabilizes the pulse. Below some 
critical value for Pco the pulses will propagate through 
both channels, essentially independent of the selected ge- 
ometry. 
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FIG. 3: (Color online) Critical angle Oc for a pulse to turn 
around a corner, as a function of dimensionless channel width 
W. For fixed CO pressure, pulses in a geometry with and W 
chosen "above" the solid lines shown are able to turn around 
the corner. The solid lines, fitted to data points, are in- 
cluded to guide the eye. The dashed and dotted lines are 
given by rc — si^^-e ) different Tc. The unit length 
of W corresponds to a real length of 3.7/im. T = 535. 5i^, 
P02 = 2.0 X 10~^mbar. 



B. Quantifying the geometry-induced instability in 
reactive pulse propagation 

We now begin a quantitative exploration of the 
geometry-induced instability of pulse propagation in our 
model system. Junctions between different linear "cor- 
ridors" are an important building block for complex ge- 
ometries. The two geometrical parameters characterizing 
a junction are the width of the channel W and the angle 
at the corner (Fig. I^J. When the reaction conditions 
are fixed, depending on the choice of W and 6 , the pulse 
may or may not be able to pass the corner, as shown in 
Fig. El 

When a pulse attempts to turn around a corner, if W 
is small, the entire pulse becomes curved as a circular 
arc due to the influence of the no-flux boundaries, and 
the curvature of the pulse is given by ^ = £l!i(^i^. if 
this geometry-determined curvature becomes larger than 
a critical value ^ , a constant determined by the dynam- 
ics of the system and the reaction conditions, the pulse 
becomes unstable and disappears j3l|. This disappear- 
ance occurs dynamically through a "decollation" mecha- 
nism 32]. When W is large, the pulse only curves locally, 
close to the corner, and further increase in W has little 
influence on this high local curvature of the pulse. When 
W is above some critical value (about 6 for curves in 
Fig. this local curvature of the pulse reaches a min- 
imum and becomes independent of W and only a func- 
tion of ^. If this minimum value of local curvature is still 
larger than ^, the pulse will decollate from the corner 
and fail to pass through. We can thus rationalize the 
existence, for each given set of reaction conditions, of a 
critical angle below which the pulse can not propagate 



around the corner for any W. 

Another geometric structure with features similar to 
the Y-j unction is the rhomb constriction, which arises in 
microcomposite Ti02/Pt checkerboards (Fig.^. 

The critical angles Oc (for the Y-junction) and (j)c (for 
the rhomb constriction) that allow the pulse to propa- 
gate through the structure are plotted as a function of 
the width VF, as in Fig. [5J under the same reaction con- 
ditions. The y-axis in Fig. for the rhomb constric- 
tion structure is chosen to be ((/)/2 + 7r/2)/7r, so that the 
two critical curves correspond (they are almost identical) 
when W is small (< 1). The non-monotonic decrease of 
(f)c for the rhomb constriction structure can be rational- 
ized by a simple geometric observation: as W increases, 
the front and back of the pulse are approaching closer 
and closer to each other at the corner (see the insets in 
Fig. [5|); this can lead to a different mechanism for the 
decollation of the pulse. 



C. Pulse manipulation with local laser heating 

1. Using the laser to assist pulse propagation 

The local increase of temperature caused by a short 
laser "burst" accelerates the local desorption of CO and 
may thus assist the propagation of an oxygen pulse. We 
can therefore use local laser heating to assist pulse prop- 
agation around corners that would be too sharp for the 
pulse to go through at isothermal conditions. 

In Fig. ini a local increase of temperature (see the high 
temperature field close to corner in Fig. Efb')) reignites 




FIG. 4: (a) A snapshot showing a microdesigned Ti02/Pt 
checkerboard. The inactive Ti02 rhombs appear black while 
pulses (little arcs in the first snapshot) are propagating on 
the light-gray Pt(llO) surface, (b) A blow-up of the checker- 
board structure. To its right, a geometry (symmetric around 
the centerline) we used in the simulations of the rhomb con- 
striction is depicted. 
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FIG. 5: (Color online) Critical angle for pulse turning around 
the corner as a function of channel width W for the rhomb 
constriction and the Y-j unction structure. For the rhomb 
constriction structure, we plot (/)c/2 + 7r/2 vs W instead of (/)c 
vs W. For rhomb constriction (Y-j unction) structures with 
values of W and 0/2 + 7r/2 (0) chosen "above" the corre- 
sponding curve, pulses are able to turn around the corner. 
The solid lines, fitted to data points, are included to guide 
the eye. Snapshots of pulse propagating in the rhomb con- 
striction structure at point a (2.7, 0.88) and b (4.0, 0.88) are 
shown in the insets. T = 535. 5i^, Pco = 4.95 x 10~^mbar, 
Poo = 2.0 X 



FIG. 6: (Color online) Using local laser heating to assist 
pulse propagation around corners. Snapshots (a)-(d) show 
the coverage of oxygen; corresponding instantaneous temper- 
ature fields before, after ((a') and (c') respectively) as well 
as during the laser shot ((b')) are plotted in the second row. 
Laser heating is turned on for a total of 0.4s (from t=5.0s to 
5.4s) and centered at the lower junction corner; the maximum 
temperature increase there is ?tK. (a) t=4.5s, (b) t=5.0s, (c) 
t=5.5s, (d) t=8.0s. T = 535.5K, Pco = 4.95 x 10"^m6ar, 
Poo = 2.0 X 10~^m6ar. 



the CO-poisoned surface, and reattaches the decollated 
oxygen pulse back towards the corner (close to the corner 
in Figs-injc) andEJd)). This result has been confirmed 
experimentally as we will see below. 



2. Using the laser to prevent pulse propagation 

Laser heating can also be used to eliminate a pulse, as 
shown in Fig. [3 First, we apply local laser heating for 
a short time to locally ignite the CO-poisoned catalytic 
surface far ahead from the pulse without creating a new 
pulse (Figs.lHb) and[7|;c)). 
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FIG. 7: (Color online) Using local laser heating to prevent 
pulse propagation. Similar to Fig. |^ instantaneous temper- 
ature fields ((a')-(c')) are plotted below the corresponding 
oxygen coverage ((a) -(c)). The local laser heating is turned 
on for a total of 0.4s (from t=1.5s to 1.9s) and centered in 
the middle of the lower channel entrance with a maximum 
temperature increase of 4.5K. (a) t=1.0s, (b) t=1.5s, (c) 
t=2.0s, (d) t=7.0s, (e) t=10.0s, (f) t=12.0s. T = 
Pco = 4.95 X IQ-^mbar, Po^ = 2.0 x IQ-^mhar. 

The local surface slowly reverts to the quenched steady 
state after the laser heating is turned off (as in Fig.[7fd)). 
There exists a refractory period before this area of the 
surface can be reignited, because the dynamics of the 
Pt phase-reconstruction are slow. If a propagating pulse 
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passes through this area before the local surface recov- 
ers, the pulse may not propagate through this area and 
"evaporates" (see Figs.EJe) andEff)). 



IV. EXPERIMENTS 
A. Experimental setup 

In earlier experiments, propagating chemical waves 
had to be visualized in the presence of inert Ti boundaries 
l3^ . These measurements were all performed using a 
photoemission electron microscope (PEEM), which has a 
sufficient lateral resolution and images the changes of the 
work function due to the adsorbed species. To achieve its 
high resolution, the sample has to be at a distance of 4 
mm to the electron cathode lens, restricting severely any 
access to the surface and, in particular, preventing laser 
beam addressing. 

At first look, ellipso-microscopy for surface imaging 
(EMSI) seems to be a more appropriate choice, since it 
leaves the whole surface totally open to additional ex- 
periments [2^. But EMSI relies on differences in ad- 
sorbed layer thickness. A microstructured surface fea- 
tures boundaries, laterally limiting the chemical waves, 
that are made out of Ti or similar metals and are sev- 
eral hundred A tall. The signals from those boundaries 
would easily saturate the image, leaving no contrast to 
observe the reactive pattern formation. A different imag- 
ing method was therefore necessary, one that provided 
good contrast for the observation of the reaction dynam- 
ics, but which was preferably insensitive to the Ti mi- 
crostructures. Additionally a sufficiently large working 
distance between the imaging instrument and the sam- 
ple was necessary to allow for local addressing of surface 
activity by means of focused laser light. 

Reflection Anisotropy Microscopy (RAM) conveniently 
combines the required properties. RAM has been exten- 
sively used to study CO oxidation on platinum. The 
first setup basically consisted of a classical ellipsometric 
configuration under almost normal incidence 25]. This 
setup was later improved by utilizing the intrinsic proper- 
ties of a Foster prism and working under exactly normal 
incidence |3^]. However, the spatial resolution of RAM 
(30/im) was insufficient to image the relatively weak and 
thin reaction pulses presented in this paper. In order to 
improve the spatial resolution significantly, a new instru- 
ment was designed and built. Figure [51 shows a sketch of 
the experimental setup. 

The Pt single crystal is located in a UHV chamber. 
As a light source for RAM the 488nm line of an Ar~^ 
laser is used. The laser light is collimated and after- 
wards polarized by means of a Glan-Thompson prism. 
The dielectric film within the polarizing beam splitter 
cube reflects about 99.9% of the incoming vertically pre- 
polarized light onto the Pt sample. During reflection 
at the Pt(llO) surface the polarization may change de- 
pending on the coverage dependent surface reconstruc- 
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FIG. 8: Schematic image of the experimental setup (not to 
scale). For details see text. 

tion l^^- The light propagates back to the beam split- 
ter. The component with retained vertical polarization 
is reflected back to the light source. The other compo- 
nent which is polarized parallelly, due to interactions with 
the anisotropic Pt(llO) surface, is transmitted and used 
for imaging of the crystal. Twice the light has to pass 
a UHV window. This window is specially designed to 
minimize stress-induced birefringence. The spatial reso- 
lution of this new setup is Sjam. Thus, also thin pulses 
close to the limit of stability can now be observed with 
RAM. Through a second UHV window, the light of an 
additional Ar~^ laser (514nm) can be focused onto the Pt 
surface. The position of the laser spot can be controlled 
via two computer-controlled galvanic mirrors. In order 
to remove unwanted 514nm stray light from the RAM 
imaging path, a dielectric filter is used. 



B. Experimental results 

In the experiments, we were able to guide pulse prop- 
agation through different local laser heating protocols. 
We can assist the propagation of a pulse in the direction 
originally prohibited by the boundary curvature-induced 
instability through a short duration laser shot (Fig. 
alternatively, we can apply excessive heating that actu- 
ally destroys the pulses (Fig.^J. 

In Fig. IHl an oxygen-rich pulse is first ignited at the 
upper left channel and propagates toward the right. It is 
able to turn around corner A, but it decollates at corner 
B unless it is assisted by local laser shot. After applying 
a laser shot at the lower end of the pulse at 5, the in- 
crease of local temperature reduces the CO coverage of 
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FIG. 9: Time- lapse images (an interval of 3s between each two 
snapshots) showing a dying oxygen-rich pulse "rescued" by a 
laser shot centered at B. The arrows indicate the propagation 
direction of the pulses. T = 4:55K, P02 — 2 x 10~^ mbar, 
P^^ = 6.42 X l^~^mhar. (a) No laser shot is applied. The 
pulse decollates at B and enters the upper right channel, (b) 
A laser shot is applied for 50ms with a power of 500mVK when 
the lower end of the pulse reaches point B. The pulse does 
not decollate at and splits into two at C entering both 
channels. 

Av>%4 aU Ay? 

•If 0,25 mm it 

FIG. 10: Guidance of CO-rich pulses in a labyrinthine struc- 
ture. The circles indicate the location where the laser shot is 
applied. The laser shot duration is Is with a power of 230mVK, 
which destroys any part of the pulse within the yellow circle. 
T = 451i^,Po2 = 2 X IQ-'^mbar, Pco = 5.50 x IQ-^mbar. 
(a) No laser shot is applied and the pulse enters all the three 
channels, (b) A laser shot is applied when the lower end of 
the pulse is at corner A of the crossing. The pulse only enters 
the upper channel, (c) A laser shot is applied when the upper 
end of the pulse passes corner B of the crossing. The pulse is 
thus prevented from entering the channel on the right. 

the local CO-poisoned surface, which makes the adsorp- 
tion of oxygen easier, so that the pulse is able to pass 
around the B corner (Fig.IHfb)). 

If excessive heating is applied to a propagating CO-rich 
pulse, removing enough CO from the local surface, the 
pulse can be locally destroyed. In Fig. [TUT b). we apply 



the laser heating for Is when the lower end of the pulse 
reaches corner A of the crossing. This erases the lower 
part of the pulse within the yellow circle and prevents 
the pulse from entering the right and lower channels. In 
Fig. irnT c). we wait until the upper end of the pulse passes 
the B corner of the crossing before applying the laser 
heating. Only the pulse propagation in the right channel 
is thus obstructed. 



V. SUMMARY AND CONCLUSIONS 

Spontaneous pattern formation, and the dynamics of 
the resulting patterns in reacting systems, can be con- 
trolled at several levels. Designing the geometry of 
the catalytic domain, and using a galvanometer mirror- 
manipulated laser beam, we have demonstrated here 
through experiment and computation, that the propa- 
gation of pulses in complex geometries can be guided, 
facilitated or forbidden in real time. Sudden and intense 
boundary curvature changes can lead to a fundamental 
"decollation" instability; on the other hand the laser- 
induced heating can either assist or prohibit pulse prop- 
agation depending on its intensity and location. 

The combination of spatially fine-grained sensing with 
spatially fine-grained actuation opens up a wide array of 
possibilities in manipulating physical processes in com- 
plex microgeometries. Here we studied heterogeneous 
reacting systems; similar tools can be used to imple- 
ment spatiotemporal networks in homogeneous reacting 
excitable media (e.g. Ref. 36]) or electrochemically re- 
acting systems (e.g. Ref. |37J). Beyond chemically react- 
ing systems, such tools are becoming increasingly useful 
-and used- in applications as diverse as droplet forma- 
tion and mixing in microfiuidics 38], directing cell mi- 
gration 39] or manipulating coherent matter- waves 4fl.]- 
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